Conditional squeezing of an atomic alignment 
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We investigate the possibility to perform Quantum Non Demolition measurements of the collective 
alignment of an atomic ensemble in the case of a -F > 1 spin. We compare the case of purely vectorial 
and purely tensorial Hamiltonians and show how to achieve conditional squeezing or entanglement 
of atomic alignment components. 
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I. INTRODUCTION 

The reduction in the quantum fluctuations of an 
atomic ensemble angular momentum has recently re- 
ceived much attention in connection with quantum in- 
formation, high-sensitivity frequency measurements and 
high-precision magnetometry. Such spin-squeezed atomic 
states may be obtained via non-linear interaction pro- 
cesses between an ensemble of A-atoms and cavity fields 
P, [3| or by direct mapping of a squeezed state of light 
onto the ground state atomic spin 0, 0, H Another 
approach consists in probing the atomic angular momen- 
tum in a Quantum Non Demolition (QND) manner in 
order to reduce the quantum fiuctuations of one of its 
components below the standard quantum noise 0, Q- 
The atomic squeezing is then conditioned on the QND 
measurement result and can be actively fed back to the 
atomic angular momentum using a magnetic field . So 
far, these protocols have been implemented with cesium 
atoms [13, [lH- Since the angular momentum is greater 
than 1/2 a full description of the atomic state not only 
requires to take into account the three components of the 
angular momentum, but also the higher order tensorial 
components. This means that one has to add to the sim- 
plified effective QND Hamiltonian three of 
the five components of the atomic alignment, which in 
general will perturb the measurement of the orientation 
Fz. It is however possible to choose the atomic detun- 
ing with the excited states such that their contribution 
is zero or negligible [13, • 

The goal of the present paper is to investigate high- 
angular momentum atom situations in which the Hamil- 
tonian is not purely vectorial and show how it is actu- 
ally possible to reaHze QND measurements of the atomic 
alignment components. Such measurements may then al- 
low for squeezing not only the quantum fiuctuations of 
the atomic orientation, but also those of the alignment, 
which are involved in several atom/light quantum inter- 
face protocols [3, For instance, by achieving condi- 
tional squeezing of the alignment of an atomic ensemble 
combined by single atomic excitation retrievals using the 
"DLCZ" protocol [ll,[lB|, it is possible in principle to pro- 
duce exotic atomic states with a non-Gaussian Wigner 



functions, in a way similar to non-Gaussian optical states 
In addition to being a tool for atomic quantum 
noise studies, controUing the fiuctuations of the atomic 
alignment may be of interest for improving the precision 
of magnetometers [lol . fiol . [20| . In order to draw simple 
conclusions we shall limit ourselves to a first order Hnear 
atom-field interaction in the optical pumping regime, but 
we note that interesting possibilities may also be offered 
by orientation/alignment conversion [2l|| and non-linear 
selective addressing of high-rank atomic polarization mo- 
ments [H. 

In Sec. [n] we give the effective Hamiltonian and derive 
the atom-field evolution equations. After reviewing in 
Sec. Uni the well-known vectorial Hamiltonian situation 
leading to QND squeezing of the orientation, we exam- 
ine in Sec. lIVI the purely tensorial Hamiltonian situation. 
We highlight the differences with the vectorial situation 
and show how QND measurements of the alignment can 
be performed, leading to conditional squeezing or entan- 
glement of the atomic components. The effect of spon- 
taneous emission losses on the obtainable squeezing and 
the experimental feasibility are discussed in Sec. [V] in the 
case of rubidium atoms. 



II. HAMILTONIAN AND EVOLUTION OF THE 
SYSTEM 



We consider an optical field propagating along z which 
interacts with an iV-atom ensemble in the low saturation 
regime and consider slow processes as compared to the 
evolution of the excited state populations and optical co- 
herences, which can be adiabatically eliminated. In this 
case, the effective Hamiltonian describing the atom-light 
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interaction can be written as [13, \2S 
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F (resp. F') is the total angular momentum of the ground 
state (resp. of one of the excited states), and its cartesian 
components are denoted by F^^y^z- fF' is the resonant 
cross-section of the F ^ F' transition, and Ai?/ is the 
probe one-photon detuning with respect to this transition 
(> if blue detuned). A is the field cross-section and L 
the length of the A^-atom medium. The vectorial and 



tensorial polarisabilities are denoted by ay and and 
their exact form, given in Refs. (isl . [23| . is reminded in 
Appendix[Al The definition for the Stokes operators used 
throughout the paper is 
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where the field a wit h frequency w is defined by = 
£o{a + at) and So = \/huj/2€oAc. 

To simplify the discussion and relate it to the experi- 
mental situation which will be considered in Sec. |Vl we 



assume in the following an = 1 total ground state spin, 
but the physical conclusions would actually remain the 
same for a higher angular momentum. The irreducible 
tensor operators for F = 1 are given by [l^l 
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with F±— Fx± iFy. In this case, the Hamiltonian reads 
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The anti-hermitic terms in the Hamiltonian of Eq. ([T]) 
are due to optical pumping. For an off-resonant in- 
teraction, these anti-hermitic terms may be neglected, 
although their contribution should be considered care- 
fully when it comes to optimizing the squeezing as it 
will be shown in Sec. [Vl If these terms are neglected 
the evolution of the atomic operators is simply given by 
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We Hmited ourselves to this set of three operators, since it is a closed system under Hint and allows for conditional 
squeezing of r| -I- or T| -I- as we will show later. The terms oc ay in Eq. lfT2|) correspond to light-shifts, and 
the ones oc ar to Raman processes involving coherences between sublevels with |Ami?| = 2. Under the slowly varying 
envelope and paraxial approximations the field evolution equations read 
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The terms oc ay in Eq. (fT3|) correspond to the well known 
Faraday rotation. In the following, we will consider the 
Stokes operators before (in) and after (out) the interac- 



tion, integrated over the pulse duration T : S*"/°"* ~ 
/p dts {0/L,t), and the collective atomic operators be- 
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fore/ after the interaction A™/°"* = dzA{z,Q/T). s 
and A have been normahsed so that S™/°"* and A™/°"* 
are dimensionless. We note that the evolution equa- 
tions (|12|13p can alternatively be deduced following the 
methods of [H, ^ for the atoms and [13, [H] for the 
photons. 



III. VECTORIAL HAMILTONIAN 

For ttT = 0, the Hamiltonian ([!]) reduces to the well- 
known "QND" Hamiltonian h^^^S^F^, which allows 
for non destructively measuring via a measurement of 
the conjugate observable of S^^ as was shown in (tI. I^. fiol. 

We briefly review the principle of this conditional 
squeezing of the orientation before generalizing it to an 
alignment in the next section. 

Prior to the measurement of Fz, the atoms are pre- 
pared in a coherent spin state oriented along x, i.e. the 
atoms are pumped into an eigenstate of F^; . The values of 
the components orthogonal to the mean spin, Fj, and F^, 
are unknown a priori, and because of the commutation 
relation [F.y, F^] = iY^ = iN, their standard deviations 
satisfy AFyAF^ > N/2. When there exist no correlation 
between the transverse components, such as in a s ample 
prepared by optical pumping, AF^ = AF^ — N/2. 
The atoms are placed in zero-magnetic field. The probe 
is linearly polarized ((S) = nx). Integrating the evo- 
lution equations, one obtains the following input-output 
relations : 
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The operators have been normalized so as to have 
unity variance when they are in coherent states {x,p = 
Fj/,z/\/A^/2 and Sy,z = Sy^z/^/n). It is clear that by 
measuring the fluctuations of s°"* one acquires informa- 
tion about the fluctuations of p {Fz is measured non- 
destructively via the Faraday rotation of the probe po- 
larisation it induces). The measurement is all the more 
accurate that the vectorial coupling strength 





in , in 


^out 


= X + KySz 






pOut 


= p'" 


out 








out 









Ky = OLv 



err 

ilA 



(18) 



is large. One therefore conditionally squeezes the atomic 
orientation. The variances of the transverse components 



after the measurement-induced projection of s™* can eas- 
ily be shown to be those of a minimal spin-squeezed state 
[23,[33,[3ll 
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FIG. 1: (a) 4-level scheme leading to a vectorial effective 
Hamiltonian Hx oc s^x. (b) 3-level A scheme leading to a 
tensorial effective Hamiltonian H/^ oc s^x + Syp. 



IV. TENSORIAL HAMILTONIAN 

A. Single-pass interaction 

Another interesting situation is the opposite case of 
a purely tensorial Hamiltonian, in which ay = 0. In 
practice, the interaction involves several hyperflne ex- 
cited states F' , so that it is possible to choose the de- 
tuning such that the various vectorial contributions van- 
ish (JF'Cty -^—^ ~ 0, while the total tensorial con- 
tribution J2f' '^F'Ctj- does not. It is then possible 
to realize a conditional measurement of the alignment in 
this particular situation. Let us assume that the atoms 
are prepared in a coherent spin state along z. The con- 
jugate transverse components in this case are '^'^'^'^-^ 



and 



smce 

malize them as previously: — (j.2 
p — (Tj — T?_2)/(iViV') and assume a circularly polar- 
ized probe: (S) — nz. 

The result of the integration of Eqs. I|12|13p can be 
found in Ref. and is reminded in Appendix [BJ It 
yields input-output relationships involving complex spa- 
tiotemporal modes for the flelds and the atoms. For a 
thin medium (kt ^ 1), they lead to the following input- 
output relations : 
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FIG. 2: Schematic of the double pass configuration proposed 
to perform a QND measurement of a collective atomic align- 
ment. 



The interaction is obviously not QND, since both com- 
ponents of the spin are now modified by the field, and 
conversely. This arises from the fact that the effective 
Hamiltonian in this case, Ha cx SxX + Syp, is quite dif- 
ferent of the previous vectorial situation Hx oc s^x, and 
now involves both quadratures (Fig. [1]). As noted in Q, 
this tensorial Hamiltonian corresponds to a linear cou- 
pling between two harmonic oscillators which, when res- 
onant, allows for efficient quantum state transfer between 
atomic and light variables and may be used in quantum 
memory protocols. As the coupling strength kt is in- 
creased, and 5°"* (and p°"* and s°"*, respectively) 
coherently exchange their fiuctuations, and it can indeed 
be shown that, when the collective coupling strength kt 
is large, the field fiuctuations are efficiently mapped onto 
the atoms and vice- versa [1, @| . 

However, since the atomic variables evolve during a 
single-pass "tensorial" interaction, it is a priori not well- 
suited for QND measurements. Nevertheless, it is still 
possible to perform a conditional measurement of the 
alignment by using two ensembles a and b (or by mak- 
ing two successive passes in one ensemble) with opposite 
mean orientations, such that Hint oc {xa + xi,)sx + {pa + 
Pb)sy and [xa + Xb, Pa + Pb] = 0. As will be detailed in 
the next Sections, this restores the QND character of the 
interaction. The physical interpretation is that both field 
quadratures are written onto the atoms in each ensemble, 
but, because of the opposite orientations, their contribu- 
tions cancel out, leaving the total afignment components 
unchanged, while the field still carries out information 
about both atomic alignment components. 
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atomic 
orientation 
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FIG. 3: Schematic of the double ensemble configuration pro- 
posed to perform a QND measurement of a collective atomic 
alignment. 



atomic medium, so that the non linear coupling allows for 
unconditional squeezing. We also note that similar ideas 
have been proposed for quantum memories and squeezing 
generation in Ref. |34| in the case of a purely vectorial 
hamiltonian. In [3j|, the second pass is used to couple 
the second quadrature of light to the atoms. Again, the 
tensorial situation is quite different, since the Hamilto- 
nian directly couples the two quadratures of light to the 
atoms. However, due to the SxX + SyP form of the Hamil- 
tonian, in order to perform a QND measurement of the 
alignment, one has to compensate for extra precession 
terms, which can be successfully done with a double-pass. 

The double-pass interaction subsequently leads to 
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s""- 1 = 1- + o(k|). For a 



The measurement of s°"* (resp. s°"* ) projects 

out' t-^rxc^-^ ^OUt' ^ 

out' I ^OUt' 

moderate value of kt = 0.35, these variances are ~ 0.5, 
significantly smaller than the standard quantum limit. A 
rigorous derivation of the conditional variance, including 
the terms of order k'^ in the input-output relations l(20l - 
[23| . can be obtained from the exact results of I|B1HB4P 
and leads to exactly the same conditional variance. The 
measurement of x performed this way is fully QND only 
for small values of kt, and will only result in a limited 
squeezing in principle. We now turn to a situation al- 
lowing for a QND measurement of the alignment for any 
value Kt- 



B. Double-pass interaction 

We first consider the double-pass geometry depicted 
on Fig. [21 A quarter- wave plate is inserted between the 
cell and the mirror, its neutral axis being aligned along 
X. We assume that the same pulse successively propa- 
gates back and forth in the atomic ensemble, with no 
temporal overlap. This is different from the situation of 
Refs. [H, [s^l, where the pulse interacts with itself in the 



C. Double-cell interaction 

Alternatively, a single-pass interaction can be per- 
formed with two atomic cells having opposite orientations 
- as in [ll| - in order to entangle the alignment compo- 
nents of two atomic ensembles. As shown in Fig. [3b the 
light pulse propagates through two ensembles (a) and (b) 
prepared with opposite orientation (F°) = — (F^) = N, 
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so that the input-output relationships now read 



[Xa + Xb) 
{Pa+PbT 



{Xa + Xh) 

{pa+PbY"' 

S™ + KT (Pa+PbY 



S™ - KT {Xa + Xb) 



(29) 
(30) 

(31) 
(32) 



If the probe pulse duration is much longer than the 
time required to propagate through the two cells, the 
pulse interacts simultaneously with the two ensembles. 
In this experimentally accessible situation, the previous 
relations hold to any order in kt, and the measure- 
ment is perfectly QND. Similarly to the vectorial situ- 
ation, measuring s°"* squeezes the variance of Xa + Xb 
to 2/ (l + 2k^) . Note that one has [xa + Xb,Pa + Pb] = 
i {Fz a + Fzb) '2/N = 0, since the two ensembles have op- 
posite orientations. It is therefore possible to squeeze not 
only the fluctuations of Xa+Xb, but also those of Pa +Pb- 
As can be seen from Eq. (|3T1) . sending a second pulse 
and detecting s°"* instead of s°"* allows for squeezing 
Pa + Pb, leaving the alignment of the two ensembles en- 
tangled. The expected value of entanglement obtained is 

AePR = A2 {Xa + Xb) + A2 {pa + Pb) = 4/ (l 2^) < 

4. Note that the result is the same as in the vectorial 
situation of [lH, but the physical situation is rather dif- 
ferent, since the tensorial situation requires a double-pass 
for the alignment measurement to be completely QND. 



V. ATOMIC NOISE AND EXPERIMENTAL 
VALUES FOR *^RB 

A. General Hamiltonian and non-zero frequency 
noise measurements 

For a single pass and in the case of a non-zero Ky, the 
input-output relations read to flrst order in Ky, kt 
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In the double-pass geometry described in Sec. IIVB| the 
vectorial contributions cancel out and Eqs. l|25H28p are 
left unchanged, so that the alignment can still be condi- 
tionally squeezed in this scheme. In the double-cell con- 
figuration, the vectorial contributions to the field evolu- 
tion (Faraday rotation) naturally cancel out, but the vec- 
torial contributions to the atom evolution (light-shifts) 
do not. However, a z-aligned magnetic field with Larmor 
frequency fii = "^KK^^^^ '^^^ compensate for these 
light-shifts. 

Another experimentally relevant issue is the measure- 
ment of the Stokes parameters fiuctuations. Techni- 
cal noise is in general smaller than the quantum fiuc- 
tuations of light only for higher-frequency components 



(typically above 0.1-1 MHz). It is therefore important 
to consider whether the schemes proposed in Sec. IIVBI 
and Sec. IIV CI can be extended to non-zero frequency 
noise measurements. In the double-cell configuration, 
it can easily be done by means of a z-aligned magnetic 
field. The Larmor precession couples x and p, but in 
the frame rotating at 2^ (il is defined by ftfi = /i-Bz, 
where /i is the magnetic moment of the ground level and 
Bz the magnetic field value). The input-output rela- 
tions (|29H32p then remain unchanged when making the 
substitution x {x)2n = x cos(2f2t) -f p sin(2ilt), 
P {p)2n = P cos(2f2i) — x sm{2Qt), etc. It is thus 
possible to measure in a QND manner the atomic opera- 
tors X2n and p2n through their imprints on the sidebands 
components of Sx or Sy. Unfortunately, this technique 
cannot be used in the double-pass configuration, since the 
magnetic field would have to be reversed between the first 
and the second pass, which is not very realistic experi- 
mentally. However, if measurements of the Stokes op- 
erators at non-zero frequency are more easily shot-noise 
limited, it was shown in that strong spin-squeezing 
could still be obtained experimentally in a zero-magnetic 
field/frequency situation. 



B. Atomic noise considerations 

We now discuss the intrinsic limitations brought by 
spontaneous emission noise in the tensorial Hamilto- 
nian case. For the sake of simplicity we study the case 
of a 1 0,1 transitions, with atoms oriented along 
z. Using the Heisenberg-Langevin evolution equations 
and the quantum regression theorem, we obtain for a 
1^0 transition (for which ay — —1/2, ar = — 1 and 



Kt = ky^/2 = kq): 
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(37) 
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(40) 



with Ea 



A 



N' 



A 



and fx,p, fs^. 



standard vacuum noise operators with variance unity. 

For the 1^1 transition, one has ay = —3/4 and ut = 
3/2, so that Kt = —Kv\f2 = k\ and similar equations 
can be derived. Choosing the detunings such that = 
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Ki = k/2 cancels the vectorial terms finally yields the 
following input-output relationships 



- a 
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with Ea = 
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Ao 



- f _i 1_\ r _ kT [n ( 1 

and e' = — + . One retrieves beamsplitter- 
like relations for the losses, similar to those of £p 
simply describes absorption of the probe caused by spon- 
taneous emission : the probe field is damped by a factor 
yjl — Ep, and some uncorrelated vacuum noise ^/£^fs^,sy 
is consequently added, as for the propagation through 
a beamsplitter with transmission y/1 — Sp. £a describes 
the symmetrical process for the atoms : the probe, be- 
cause of spontaneous emission, induces optical pumping 
towards a z-aligned coherent spin state (which is similar 
to mixing the probe with some vacuum). A difference 
with the vectorial situation is the presence of small con- 
tamination terms a e'. They also correspond to optical 
pumping processes which tend to align x and p along 
and Sy. To minimize the effect of spontaneous emission 
noise, one has to choose n ^ in order to have Ea ^ £p, 
as in Ref. Finally, the total spontaneous emission 
contribution in the double-pass or double-cell configura- 
tions is finally obtained by doubling Eq, Ep and e' in the 
above equations. 



C. Experimental values for Rb 

Based on these considerations we discuss the values of 
squeezing or entanglement that can be expected in ex- 
periments with ^''Rb. We assume an interaction on the 
D2 fine with the atoms in the F = 1 ground state. For 
room temperature vapor cells, taking into account the 
Doppler broadening, no detuning allows for completely 
canceHng Ky. On the contrary, for cold atoms with negli- 
gible Doppler broadening, Ky can be canceled for a probe 
laser blue-detuned by Aq = 38 MHz from the F' = ex- 
cited level (red-detuned by 34 MHz from F' = 1 and 191 
MHz from F' = 2). For typical values for the density 
and volume (10^^ cm~'^ and 0.5 mm'^) of a cold atom 
cloud produced using a magneto-optical trap (the lat- 
ter being switched off during the measurement), lead- 
ing to = 0.5 X 10^ and taking a pulse of intensity 1 
/iW and duration 0.5 /iS containing n ~ 0.5 x 10^ pho- 
tons (saturation parameter ^ 10^"^, considering a cross- 
section A = 1 mm^), the previous calculations predict 
kt ^ —0.42 and ^ ~bdB of squeezing in the quantum 
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FIG. 4: Vectorial and tensorial coupling strength k,v (dot- 
ted) and Kt (dashed), and amplitude of the noise added to 
the probe Sp (plain) as functions of the normalised detuning 
= Ao/(r/2) between the probe and the -F = 1 ^ F' = 
transition of *^Rb D2 line. The experimental parameters are 
detailed in Sec lVC] The insert zooms on the detuning area 
where the Hamiltonian is purely vectorial. 



for longer probe pulses, provided that the duration of 
the pulses remains smaller than the relaxation time of 
the Zeeman coherence, or by the use of a dipole trap to 
increase the optical depth [lj|. For these parameters, in a 
double-pass or double-cell configuration, Sa = Sp — 0.14 
and e' < lO^'' (the contribution of the F' = 2 level is 
~ 30 times smaller than those of the F' = 0,1 levels, 
and is not considered). As the fiuctuations are predicted 
here to be reduced by a factor smaller than ~ 1/0.14, the 
noise added by spontaneous emission can be neglected. 

For the sake of comparison, we now discuss the relative 
strengths of tensorial and vectorial conditional measure- 
ments. In atomic vapors close to room temperatures, 
the detuning is usually chosen bigger than the doppler 
broadening in order to avoid absorption [ll|. It implies 
that the detuning has to be large as compared to the hy- 
perfine structure, and, since one has ap'O-T = for 
alkali atoms, it means that kt oc crpiaT /A^ ^ 0, 
i.e. the effective Hamiltonian is then almost purely vec- 
torial. This situation is obviously much more favorable 
for orientation than for alignment squeezing. 

However, for a doppler-free medium, it is possible to 
reduce the detuning while maintaining a small absorp- 
tion. In this case, as can be seen from Fig. IH both kt 
and Kv (and hence alignment and orientation squeezing) 
may have similar values. To compare these values, we 
consider the case of a purely tensorial Hamiltonian (i.e 
KV = 0, obtained for Aq = 38 MHz, A^^ = 13.2), and the 
case of a purely vectorial one (i.e kt — 0, obtained for 
Ao = 222 MHz, A^ = 77). In the first one, for the ex- 
perimental parameters given above, kt = —0.42, kv = 
and Ep = Ea = 0.14 whereas in the second , kt = , 
Kv = 0.03 and Ep = Ea = 0.01. This shows that the 
common idea that the vectorial coupling strength is big- 
ger than the tensorial one is not necessarily true for cold 
atom samples when the hyperfine structure is taken in 
account. 
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VI. CONCLUSION 

We have shown how to perform a QND measurement 
of a collective atomic ahgnment. This extends the pos- 
sibility to manipulate high-angular momentum compo- 
nents of a collective spin beyond the vectorial Hamilto- 
nian interaction commonly used so far in experiments 
0, [13, [lH ■ Noticeable physical differences are found be- 
tween the purely vectorial Hamiltonian situation and the 
tensorial situation. In particular, if it had been noted 
in previous work 0, @| that the tensorial situation may 
lead to coherent atom-field quantum state transfer and 
storage, we have shown here that it also allows for per- 
forming a QND measurement of the atomic alignment, 
provided that two ensembles or two successive passes are 
used. Substantial conditional squeezing values are still 
predicted for realistic experimental situations with cold 



atomic samples. We also note that these measurements 
can be used to continuously control the atomic spin fluc- 
tuations via feedback The different feedback mech- 
anisms that may be used to squeeze an atomic alignment 
will be presented elsewhere. 
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APPENDIX A: POLARIZABILITY 

The polarisabilities ay and are given by 



p, _ 3(2J^ + 1) f 2F-l ^p, 2F + 1 ^p, 2F + 3 ^p, 

"V - 2(2F'-H1)(2J+1) F ^^-^ FiF+iy^ ^ F+l^^+' 

p. _ 3(f +1)(2J- + 1) (l^p, 2F+\ ^p, , 1 p, 

" 2(2^^' + l)(2J+l) U F{F + 1)^ ^ F + l^+' 

where 5p, is Kronecker's symbol. The resonant cross-section between two levels with an isotropically populated ground 
state is 

2(2J+l)(2i^' + l) r ./' 1 J 

Cfp' — '^2level ^ { F I F' 

with (72ievei = ■ The commutators between the irreducible tensorial operators are [23 |: 

fci fca K 



F F F 
^3 



K (^s)' ^^.^ iPo)] = E i-^)""^^^' v/(2fci + l)(2fc2 + l) 

K,Q 

( fcifc2<zi<?2,ifg) (1 - (-i)^^+^^+^)T^(^^g) 



APPENDIX B: TENSORIAL SITUATION : of Eqs. (|12m3p yields [1] 

SOLUTIONS OF THE EVOLUTION 
EQUATIONS ()12m3|l 



KT [ dz f dz'x'"{z') 
Jo Jo 



J\[2kt\/z — z') 



(Bl) 



\\iz r dz'p-iz')^i^2^^f^^ 

Jo Jo Vz-z' 



We assume a single-pass interaction with ay = 0, as i , ^ 

in Sec. |IVA[ After changing the spatiotemporal frame -^^ / dts™(t) ( / dz Jo(2kt ^^(1 -t)) 
(Zjt) ^ (z = z,t = i — z/c) and making the system di- Jo \Jo 

mensionless [z^t) — > (z = z/L,t = t/T), the integration (B2) 
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and symmetrical equations for the fields 

= - KT dt dt (t )- 



JO 



^out An 



(B3) 

, Ji(2KT\/t - t') 



-KT^ d2.x'"{z) (^J^ dtJo(2KTv/t(l -z)) 



(B4) 



where Jo and Ji are the standard first order Bessel func- 
tions and the operators have been normalized so as to 
have unity variances when in coherent states. At first 
order in kt, one retrieves Eqs. 
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